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I.  INTRODUCTION 


Atmospheric  turbulence  effects  are  severe  for  many  laser  systems 
operating  in  the  atmosphere.  Unless  some  sort  of  adaptive  optics  corrections 
are  employed,  these  effects  will  greatly  reduce  the  viability  of  such  systems. 
In  conventional  phase-conjugate  adaptive  optics,  a  beacon  signal  that  origi¬ 
nates  at  the  object  of  the  imaging  system  or  at  the  aim  point  of  the  laser 
transmitter  system  is  used  to  provide  information  on  the  atmospherically 
induced  phase  errors  associated  with  the  desired  propagation  path.  This 
instantaneous  phase  of  the  beacon  signal  measured  by  appropriate  wavefront 
sensors  is  used  to  make  the  adaptive  optics  correction  by  conjugating  this 
phase . 

In  a  previous  paper,  Yura  and  Tavis^  discussed  in  some  detail  the  serious 

consequence  of  using  a  full  aperture  centroid  tracking  (i.e.,  Hartmann-type 

sensor)  to  measure  overall  wavefront  tilt  when  another  subsystem  (figure 

sensor)  measures  all  higher  order  phase  distortions.  The  figure  sensor  and 

deformable  mirror  are  used  to  remove  all  but  tilt  defined  in  terms  of  the 

phase,  weighted  over  the  aperture.  The  centroid  tracker  and  steering  mirror 

remove  tilt  as  defined  by  the  gradient  of  phase,  averaged  over  the  aperture. 

The  remaining  residual  tilt  is  what  leads  to  centroid  anisoplanatism .  This 

residual  tilt  leads  to  a  large  reduction  of  energy  on  target,  namely  Strehl 

ratios  of  the  order  of  0.03  for  D/r  =  100. 

o 

In  Ref.  1,  Yura  and  Tavis  also  discussed  the  effect  of  scintillation  on 
the  centroid  anisoplanatic  error  but  did  not  provide  numeric  estimates.  The 
purpose  of  this  report  is  to  provide  numerical  estimates  of  the  effects  of 
scintillation.  It  is  shown  that  these  are  effects  of  secondary  importance  to 
centroid  anisoplanatism. 
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II.  SCINTILLATION  EFFECTS  ON  CENTROID  ANISOPLANATISM 


We  assume  that  the  log-amplitude  fluctuation  x  is  not  zero,  that  it  is 

2 

normally  distributed  with  a  variance  a  <  1,  and  that  the  Rytov  approximation 
9  X 

is  valid. “ 

2 

In  Ref.  1  it  was  shown  that  the  tilt  variance  <0^>  and  the  tilt-centroid 
correlation  were  independent  of  scintillation  effects.  Further  it  was  shown 
that  the  centroid  variance  could  be  written  as 


<e2  (x )>  =  <9 2>  S 

c  c 


(1) 
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where  <9  >  is  the  mean  square  centroid  angle  variance  given  by 


<0  2> 
c 


<02  > 

co 


(l-€2) 


[a  +e11/3> 


»  2,r17wl  K 

2  e 


r (8/3) 


jFjCp  -  5/6,  2,  e2H  (2) 


<02  >  =  13.39  (D/r  )5/3(kD)”2, 
CO  o 


(3: 


D  is  the  optics  diameter,  rQ  the  lateral  coherence  length,  and  k  is  the  wave 
number  (=  2tt/X  where  X  is  the  wavelength).  The  quantity  e  is  the  obscuration 

3 

ratio,  T  is  the  gamma  function,  and  2^1  is  the  hypergeometric  function  . 

In  Eq .  (1),  S  is  given  by 

/  dr  rK(r)[v2D(r)  -  4V2B  (r)  -8(vB  (  r ))  2|exp[  4B  (r)] 

s  .2 - 5 - : -  O 

/  dr  rK( r )  V  D(r) 

0 

where 


D(r) 


6.88 


(r/r  ) 
o 
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(5' 
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---  -  =  1  d  R  w(R  +  0. 5r) W(R  -  0.5r) 

|  /  d2R  W(R)  |2 

+ 

W(r)  is  the  aperture  function.  For  e  =  0  and  a  uniform  aperture 


K(x)  =  [cos  *x  -  x(l-x2)  ], 


The  log-amplitude  and  the  log-amplitude  phase  correlation  functions 
may  be  written  in  terras  of  the  corresponding  structure  functions The 
explicit  functional  dependence  of  terms  in  Eq.  (4)  are  given  below. 


We  have^ 


Vr>  ’  <*2>  -  7  D/f> 


V>  -  <*x> 


D  (r)  =  0-5  (d]l (r)  +  Re  D2(r)  ) 


=  -  !m  D2(r) 


where 


DX(P)  =  ClPJ 


0.033  w2  l2  r(I)  =0 

"  I  <«)■>* 


02(P)  =  c2  /  cn  (*>  x3/ftdx[1F1(-5/6,  1,  , 


0  .033 TT 


7/6  -5/6 


.  2V  ~  2  0.041  k  ' 

<*  11 - n7F— 

cos  9 


/  C2  (x)x5/6dx 
0  n 


The  profile  is  plotted  in  Fig.  3,  Ref.  4,  and  ^F^  is  the  confluent 
hypergeometric  function.  We  also  note  for  convenience  that  the  wave  structure 

function  Eq.  (5)  may  also  be  written  as  [i.e.,  Eq.  (11)].  To  calculate  the 

integrals  appearing  in  Eq .  (4),  the  gradient  ( $)  and  Laplacian  (V2)  of  and 
D2  are  needed.  These  expressions  are  given  in  Appendix  A.  To  proceed  further 
with  the  calculation,  we  use  the  integral  expressions  for  the  confluent 
hypergeometric  functions  for  the  magnitude  of  the  argument  less  than  30  and 

its  asymptotic  expansion  for  values  of  the  argument  greater  than  30.  These 

expressions  are  given  in  Appendix  B.  The  asymptotic  expressions  are  valid 
over  most  of  the  integration  range,  indicated  in  Eqs.  (4)  to  (12).  For  high 
altitudes  and  small  values  of  r  [Eq.  (4)],  the  integral  expression  (B-5)  was 
used  to  calculate  the  functional  values  at  1000  points,  and  a  linear  extrapo¬ 
lation  was  used  between  specified  arguments.  The  expression  S  has  been 
calculated,  and  the  quantity  S-l  is  given  in  Table  1  as  a  function  of 
obscuration  ratio. 

Table  1.  Scintillation  Effect  on  Centroid  Variance 


c  (obscuration  ratio) 


STREKL  RATIO 


The  correction  factor  to  the  Strehl  ratio  may  be  found  by  replacing  the 
first  factor  on  the  right-hand  side  of  Eqs.  (34)  and  (41),  Ref.  1,  by 

<6  2>S 

[i  +-T-]  * 

<9te> 

where  S  is  given  by  Eq.  (4)  above.  When  this  is  done,  the  correction  to  1.8f(e) 
(Table  2,  Ref.  1)  is  in  the  4th  or  5th  decimal  place  for  the  large  D/r Q  (~100) 
used  in  the  calculations. 
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III.  CONCLUSIONS 


We  have  outlined  the  calculation  of  the  effect  of  scintillation  on 
centroid  an i sopl anat i sm.  We  have  found  that  the  scintillation  effect  is  three 
orders  of  magnitude  smaller  than  the  centroid  an i soplanat i c  effect  and  justi¬ 
fies  ignoring  its  contribution.  This  does  not  mean  that  scintillation  can  be 
completely  ignored  in  the  measuring  of  the  phase  for  the  adaptive  optics 
correct  ion. 
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APPENDIX  A 


APPENDIX  B 


Integral  and  asymptotic  representations  of  the  confluent  hypergeometric 


function  are: 


r(b-a)  T(a) 
r(b) 


a-1,,  .b-a-1 


F  (a,b,z)  =  /  exp(zt)  t  (1-t) 

0 


(B-l) 


Re  b  >  Re  a  >0 


^F^(a,b,z)  ^  exp(+iwa)z 

f(b)  =  r(b-a) 


R-l  (a  ) ( 1+a-b) 

{  y  _2 —  — 

n=0  n- 


n  (-z)‘n+  o(|z!K)} 


,  .a-b  s-1  (b-a)  (1-a)  „ 

exp(z)  {  j  - » - a  (*)-"  +  0(|z!R  )) 

;  n=0 


(B-2) 


aQ=l,  a1=a,  an=(a+n-l ) •  •  *(a) 


1  3 

|  z  |  large:  the  positive  sign  if  -  y  ir  <arg  z  <^-  it 

3  1 

the  negative  sign  if  -  ^  tt  <arg  z  tt. 


Note  that  the  restriction  on  Eq.  (B-l)  means  that  ^F^(-  -g-  ,  1,  z)  does  not 
have  an  immediate  integral  representation.  Instead,  the  recursion  relation  is 
used : 


bM(a,b,z)  -  bM(a-l,b,z)  -  zM(a,  b+l,z)  =  0 


(B-3) 


to  find 


lFr“  6’  lf  lMh*  l*  ~  z  l  Fl  (fj  ’  2‘ 


(B-4) 


The  needed  integral  representations  are  given  by 


1.  z)  =  — 7  1  s  ■  /  exp(zt)t  ^(1-t)  1/,&dt  , 

hi)r(f)  o 


lFi(r*  2»  z )  = — r  -1-- f  exp(  zt)  c  5/^(l  -t)  5//6dt  , 


iFi(r»  3»  z)  = — — —  f  exp(zt)t 1  ^(1  -t)  5^dt  . 

6  o 


The  highest  order  terras  in  the  asymptotic  expansion  are  given  by 


f  5  .  5/6 

„  r  5  ,  ,  ~  exp^ 

1F1^  6’  l»  Z>  t(  11/6] 


[l  -  (5/6)2z_1  +  0.5*(1/6)2(5/6)2z  2] 


1  ~  exp(r  iri)z  6  as  t 

i^i*  2*  z]  "  — tttt: — ^  +if  2  2~2i 

^  6j 


+  exp( z)z 

T h|) 


-1  1/16 


iFi(f  3>  :  2  •  I 


r  7  .s  -7/6 

exp[7,l)z  .  15  -1  .  ,.4  55  -2, 
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